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HE: [H 49] SEXE EA Ascosphaera apis ( 简称 球 淖 菌 ) 是 一 种 特异 性 侵 染 蜜蜂 幼虫 肠 道 的 致死 性 
真菌 病原 。 微 小 RNA(microRNA, miRNA) 是 一 种 在 转录 后 水 平 对 mRNA 进行 负 调控 的 关键 调控 
因子 。 本 研究 旨 在 对 意大利 蜜蜂 Apis mellifera ligustica ( 简称 意 蜂 ) 幼虫 肠 道 在 球 吉 菌 胁迫 后 期 的 
差异 表达 miRNAsS(DEmiRNAs) 及 其 靶 基 因 进 行 深 入 分 析 , 为 揭示 DEmiRNAs 在 胁迫 应 答 中 的 作用 
提供 重要 信息 。【 方 法 】 利 用 small RNA-seq(sRNA-seq) 技 术 对 正常 及 球面 菌 侵 染 的 意 蜂 6 Hg db 
忠 肠 道 ( 分 别 表 示 为 AmCK 和 AmT) 进 行 测序 ,通过 相关 生物 信息 学 软件 对 DEmiRNAs 进行 预测 和 
分 析 。 利 用 TargetFinder 软件 预测 DEmiRNAs 的 靶 基 因 ,然后 利用 Blast 软件 对 靶 基 因 进 行 CO 和 
KEGG 数据 库 的 功能 注释 。 通 过 Cytoscape 软件 构建 DEmiRNAs 与 其 靶 mRNAs 的 调控 网 络 。 通 过 
茎 环 实时 荧光 定量 PCR(stem-loop RT-qPCR ) 验证 测序 数据 的 可 靠 性 。【 结 果 】AmCK 和 AmT 的 测 
序 分 别 得 到 9 230 496 和 10 823 667 条 有 效 序列 标签 ;AmCK vs AmT 比较 组 中 包含 15 个 上 调和 6 
个 下 调 miRNAs, 分 别 结合 3 503 和 3 252 个 靶 基 因 ,它们 可 分 别 注释 到 40 和 39 个 GO terms 以 及 
104 和 99 条 代谢 通路 ;进一步 分 析 发 现 调控 网 络 中 的 17 个 DEmiRNAs 靶 向 结合 116 个 与 丝氨酸 蛋 
白 酶 相关 的 mRNAs ,14 个 DEmiRNAs 3&5] 25 4- 54 个 与 泛 素 介 导 的 蛋白 水 解 相 关 的 mRNAs, Stem- 
loop RT-qPCR 验证 结果 显示 随机 选择 的 4 个 DEmiRNAs (miR-251-x，miR-9277-y，miR-1672-x 和 
miR-4968-y) 的 表达 量变 化 趋势 与 测序 结果 一 致 ,表明 测序 数据 真实 可 信 。[【 结论 ] 本 研究 率先 对 意 
蜂 幼 忠 肠 道 在 球 淖 菌 胁 连 后 期 的 DEmiRNAs 及 其 靶 基 因 进 行 预测 与 分 析 , 提 供 了 miRNAs 的 表达 
谱 和 差异 表达 信息 。ame-miR-927b，miR-429-y 和 miR-8440-y 等 DEmiRNAs 可 能 参与 对 丝氨酸 蛋 
白 酶 的 调控 和 对 泛 素 介 导 的 蛋白 水 解 的 调控 ,DEmiRNAs 与 靶 基 因 之 间 存 在 复杂 的 调控 网 络 关系 。 


DEmiRNAs 可 参与 到 意 蜂 幼虫 肠 道 与 球 喜 菌 间 的 互 作 。 
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Abstract: [ ^im] Ascosphaera apis is a deadly fungal pathogen that specially infects the larval gut of 
honeybees. MicroRNA ( miRNA) plays a key negative role in regulation of mRNA at the post- 
transcriptional level. This study aims to comprehensively analyze the differentially expressed miRNAs 
( DEmiRNAs) and their target genes in the larval gut of Apis mellifera ligustica during the late stage of A. 
so as to provide important information for revealing the roles of DEmiRNAs in the stress 
response process. [ Methods] Normal and A. 


(represented as AmCK and AmT, 
technology, followed by prediction and analysis of DEmiRNAs. 


apis stress, 
apis challenged 6-day-old larval gut of A. m. ligustica 
respectively) were sequenced using small RNA-seq ( sSRNA-seq ) 
Target genes of DEmiRNAs were 
predicted with TargetFinder, and annotations of target genes in GO and KEGG databases were performed 
using Blast. The regulation networks between DEmiRNAs and the target mRNAs were constructed using 
Cytoscape. The reliability of the sequencing data was verified by stem-loop RT-qPCR. [Results] Deep 
sequencing of AmCK and AmT samples respectively generated 9 230 496 and 10 823 667 clean tags. 
There were 15 up-regulated and 6 down-regulated miRNAs in AmCK vs AmT comparison group, binding 
3 503 and 3 252 target genes, which could be annotated to 40 and 39 GO terms as well as 104 and 99 
metabolic pathways, respectively. Further investigation indicated that 17 DEmiRNAs within the regulation 
networks could bind to 116 serine protease-related mRNAs, and 14 DEmiRNAs were capable of binding 
The result of stem-loop RT-qPCR showed that 
miR- 
indicating the 


to 54 ubiquitin-mediated proteolysis associated mRNAs. 
the variation trends in the expression levels of randomly selected four DEmiRNAs ( miR-251-x, 
9277-y, miR-1672-x and miR-4968-y) were consistent with those in sRNA-seq data, 
credibility and reliability of our sequencing data. [Conclusion] Our study takes the lead in the prediction 
and investigation of DEmiRNAs and corresponding target genes in the larval gut of A. m. ligustica during 
the late stage of A. apis stress, and provides expression patterns and differential expression information of 
miRNAs. The DEmiRNAs including ame-miR-927b , miR-429-y and miR-8440-y may be involved in the 
regulation of serine protease and ubiquitin-mediated proteolysis. Complex regulations exist between 
DEmiRNAs and their target genes. DEmiRNAs are believed to participate in the interactions between the 
larval gut of A. m. ligustica and A. apis. 

immune defense; 


Key words: Apis mellifera ligustica; Ascosphaera apis; stress response; microRNA; 


target gene 


2E EPRE A Ascosphaera apis( HEREKE VN) 是 一 








微小 RNA (microRNA, miRNA) 是 一 类 长 度 约 





种 可 导致 蜜蜂 幼虫 震 患 白垩 病 的 特异 性 真菌 病原 。 
患 病 幼 虫 的 躯体 先 肿胀 微软 ,后 失 水 缩 成 坚硬 的 块 
状 物 并 在 预 晴 期 内 死亡 ,从 而 削弱 蜂 群 群 势 ,给 各 国 
养 蜂 生产 带 来 较 大 损失 ( 赵 红 霞 等 , 2014) 。 我 国 养 
蜂 生 产 中 的 常用 蜂 种 为 意大利 蜜蜂 Apis mellifera 
ligustica ( 简称 意 蜂 ) 和 中 华 蜜蜂 Apis cerana cerana 
(简称 中 蜂 ) ,前 者 是 白垩 病 的 易 感 蜂 种 , 常 大 范围 
其 发 白垩 病 。 在 协同 进化 过 程 中 ,以 蜜蜂 为 代表 的 
昆虫 进化 出 了 以 细胞 免疫 和 体液 免疫 为 核心 的 多 种 
防御 途径 来 抵御 外 来 病原 体 的 侵 染 ,其 中 蜜蜂 的 细 
胞 免疫 包括 吞噬 作用 、 黑 化 作用 等 ,体液 免疫 包括 
Jak-STAT, Toll, Imd 和 JNK 等 免疫 信和 号 通路 (Evans 
，2006)。 目 前 ,蜜蜂 对 球 圳 菌 胁迫 应 答 的 组 学 
研究 较为 滞后 ,免疫 应 答 机 制 尚 不 明确 。 
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为 18 ~25 个 核 苷 酸 的 非 编 码 RNA (non-coding 
RNA, ncRNA) ,通过 抑制 或 降解 mRNA 广泛 参与 植 
物 动物 和 微生物 转录 后 水 平 的 各 类 生物 学 调控 过 
程 (Asgari，2013; Ojha et al., 2016; Liu et al., 

2017) , miRNA 作为 关键 的 效应 因子 在 宿主 与 病原 
相互 作用 过 程 中 也 发 挥 重 要 功能 
2015)。 前 人 研究 发 现 ,在 被 水 疱 性 口 炎 病毒 
(VSV) 感 染 的 小 鼠 体 内 ,miR-24 和 miR-93 能 够 下 
调 病毒 mRNA 的 表达 水 平 (Otsuka et al., 2007) 。 
在 灰 霉 菌 Botrytis cinerea 胁迫 的 番茄 中 , sly-miR- 
482b 可 通过 靶 向 调控 4 个 植物 抗 病 蛋白 ,从 而 抑制 
外 源 病菌 的 侵 染 ( 张 伟 , 2014) miRNAs 同样 参与 
了 昆虫 与 病原 互 作 过 程 中 的 基因 表达 调控 ,如 miR- 
981 可 通过 抑制 黑 腹 果 蝇 Drosophila melanogaster 的 





























( Lucas et al., 
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Imd 信号 通路 中 抗菌 肽 分 子 的 表达 以 应 对 大 肠 杆菌 
Escherichia coli 的 侵 染 (Li et al., 2017) ;在 被 质 型 
多 角 体 病毒 BmCPV 感染 的 家 蚕 Bombyx mori 幼虫 
体内 , bmo-miR-278-3p 表达 量 显 著 下 调 , 过 表达 
bmo-miR-278-3p 可 显著 抑制 其 靶 基 因 IBP2 ,并 促进 
病毒 的 复制 (Qian et al., 2018) 。 然 而 到 目前 为 止 ， 
在 miRNA 水 平 研究 蜜蜂 与 其 病原 互 作 的 报道 信息 
极为 有 限 。 

笔者 所 在 课题 组 前 期 已 通过 二 代 测 序 技术 和 趋 
势 分 析 方 法 在 mRNA 水 平 揭示 了 意 蜂 幼虫 肠 道 响 
应 球 圳 阔 胁 迫 的 免疫 应 管 (Chen et al., 2017) ,并 初 
步 解析 了 病原 与 宿主 之 间 的 互 作 ( 陈 大 福 等 ， 
2017). HRT, Xx 95 27] E TE EK 3E DSL II GR SERE P D 
miRNAs 表达 谱 和 差异 表达 信息 仍然 缺失 , 有关 
miRNAs 在 宿主 的 免疫 防御 过 程 中 的 作用 仍 不 明 
了 。 为 在 miRNA 水 平 探究 意 蜂 幼虫 肠 道 对 球 寺 菌 
的 胁迫 应 答 ,本 研究 利用 small RNA-seq( sRNA-seq) 
技术 对 正常 及 球 洁 阔 侵 染 的 意 蜂 6 日 龄 幼虫 肠 道 进 
行 测序 ,通过 生物 信息 学 分 析 方法 对 意 蜂 幼虫 肠 道 
TE EK 3E p M id Jn] Hj HU 22 5B 46 35 miRNAs 
( DEmi RNAs) Jt FCRI AE I ETT IR A AT Br , DESEE 
de pt T SEI dd TE BR 2E R É JE: BR. mi RNAs 
AXE RA B. AH DEmiRNAs 在 胁迫 
应 答 中 的 作用 打下 了 基础 。 




















































































































1 材料 与 方法 


1.1 供 试 生物 材料 
本 研究 使 用 的 意大利 蜜蜂 幼虫 取 自 福建 农林 大 
学 蜂 学 学 院 教 学 蜂 场 ; 球 吉 菌 菌株 由 福建 农林 大 学 
蜂 学 学 院 蜜蜂 保护 实验 室 保存 与 活化 。 
1.2 球 圳 菌 培养 .活化 与 孢子 悬 液 制备 
将 实验 室 保存 的 球 圳 菌 菌株 接种 到 马 铃 蓝 葡萄 
糖 琼脂 培养 基 , 培养 亚 封 口 后 置 于 33 +0. 50 的 生 
化 培养 箱 中 培养 8 ~ 10 d 后 ,培养 基 的 中 心 区 域 可 
产生 较 多 的 黑色 孢子 吉 ,继而 将 其 刮 至 锥 形 瓶 内 ,加 
适量 无 菌 水 稀释 ,4% 冰箱 放置 10 min, 取 出 后 室温 
静 置 20 min, 加 入 干净 转子 ,使 用 磁力 搅拌 器 搅拌 ， 
直至 显微镜 检 孢 填 内 的 孢子 球 释 放 分 散 为 单个 抱 
子 , 利 用 血球 计数 板 对 孢子 进行 计数 ,纯化 的 孢子 保 
存 于 4% 冰箱 备 用 。 
1.3 ”人工 饲养 意大利 蜜蜂 幼虫 及 球 圳 菌 接种 感染 
参照 王 倩 等 (2009 ) 的 方法 , 按 蜂 王浆 63% 、 蜂 
蜜 6% .无 菌 水 30% .酵母 1% 的 质量 比例 配制 意 蜂 




































































幼虫 的 饲料 。 选 择 群 势 较 强 的 健康 意 蜂 蜂 群 ,利用 
蜂王 产 卵 控制 器 将 蜂王 内 禁 在 工蜂 清理 1 d 后 的 空 
脾 上 产 卵 ,自然 产 卵 6 h 后 将 蜂王 移出 控制 器 。 将 
集 脾 迅速 提 至 实验 室 , 结合 戴 氏 定律 ( 李 兆 英 ， 
2013) ) 和 养 蜂 生产 经 验 计算 幼虫 发 育 日 龄 ,将 2 日 龄 
幼虫 移 至 已 预 置 50 pL 饲料 的 24 孔 细胞 培养 板 中 ， 
置 于 35 +0.5%C 、 相 对 温度 (RH)90% 的 培养 箱 中 饲 
养 。 配 制 终 浓 度 为 1 x 107 个 孢子 /mL 的 饲料 , 饲 喂 
处 理 组 3 日 龄 幼虫 , 24 h 后 饲 喂 不 含 球 早 菌 的 正常 
饲料 ,对 照 组 幼虫 饲 喂 不 含 孢 子 的 正常 饲料 ,每 隔 
24 h 更 换 饲 料 。 
1.4 测序 样品 准备 及 Hlumina 测序 

球 吉 菌 孢子 被 意 蜂 3 日 龄 幼虫 摄 入 后 进入 中 
肠 , 此 时 中 肠 与 后 肠 隔 绝 ,孢子 处 于 低 水 平 萌发 状 
态 , 至 幼虫 到 预 肾 的 过 渡 期 ,中 上 肠 与 后 肠 连通 , 抱 子 
随 食 物 进入 后 肠 , 此 时 孢子 剧烈 萌发 , 菌 丝 大 量 生 
长 , 穿 透 肠 壁 后 继续 生长 ,直至 穿 透 体 壁 ,最 终 蔓延 
包 右 幼虫 全 身 ,形成 白垩 病 虫 尸 ( 陈 大 福 等 , 2017 ) 。 
考虑 到 球 圳 菌 接种 感染 前 期 (4 日 龄 和 5 日 龄 幼虫 ) ， 
孢子 在 幼虫 中 肠 仅 处 于 低 水 平 萌发 状态 ,其 与 宿主 间 
存在 互 作 ,但 并 不 剧烈 ,而 感染 球 圳 菌 的 6 日 龄 幼虫 
尚未 表现 出 明显 的 白垩 病症 状 ,但 已 临近 白垩 病 暴 
发 , 球 吉 菌 季子 在 幼虫 肠 道 内 已 进一步 萌发 ,因而 此 
时 宿主 和 病原 间 存 在 更 为 密切 的 互 作 , 窒 主 在 胁迫 
后 期 的 应 答 水 平 较 强 ,在 转录 组 水 平 能 体现 出 较 大 
差异 ,可 挖掘 出 有 价值 的 信息 和 线索 。 因 此 ,本 研究 
选择 意 蜂 6 日 龄 幼虫 肠 道 作为 测序 材料 。 

在 超 净 台中 分 别 剖 取 处 理 组 和 对 照 组 6 日 龄 幼 
虫 整个 肠 道 ,每 9 条 幼虫 肠 道 置 于 1 个 RNA-Free 的 
EP 管 中 , 液 握 速冻 后 存放 在 -80Y 的 超低温 冰箱 中 
待 用 。 本 实验 分 别 进行 3 次 生物 学 重复 。 处 理 组 编 
号 分 别 为 AmT-1，AmT-2 和 AmT-3 ,对 照 组 编号 分 
别 为 AmCK-1, AmCK2 和 AmCK-3。 委 托 广 州 基 迪 
奥 生物 科技 有 限 公 司 对 上 述 肠 道 样品 进行 单 端 测 
序 , 测 序 平台 为 Humina MiSed。 原 始 数据 已 上 传 
NCBI SRA 数据 库 , BioProject 号 : PRJNA408312 。 
1.5 测序 数据 的 质 控 及 DEmiRNAs 的 预测 

参照 本 课题 组 前 期 建立 的 方法 ( 郭 寄 等 , 2018) 
对 测序 数据 进行 质 控 和 评估 ,得 到 过 滤 后 的 miRNAs 
的 非 注释 标签 序列 (unannotated tags ) ,并 将 其 通过 
Bowit 软件 向 西方 蜜蜂 参考 基因 组 (assembly Amel 
4.5) 进 行 映射 (mapping) ,得 到 相关 tags 在 参考 基因 
组 上 的 位 置信 息 , 即 为 mapped tags。 继 而 将 mapped 
tags 通过 miRDeep2 软件 与 miRBase 数据 库 中 的 





























































































































52 昆虫 学 报 Acta Entomologica. Sinica 62 35 





miRNAs 前 体 序列 进行 比 对 ,得 到 可 能 的 前 体 序列 同 
时 鉴定 miRNAs 的 表达 ,然后 通过 TPM (tags per 
million ) 算 法 公式 (TPM =T x 105/N,T 代表 单个 miRNA 
的 tags, N 表 代 表 总 miRNAs 的 tags) 对 所 有 miRNAs 
的 表达 量 进行 归 一 化 处 理 。 各 样品 之 间 的 相关 性 系 
数 利 用 R 软件 计算 。DEmiRNAs( AmCK vs AmT) 的 
筛选 标准 为 llog, Fold changel z1 H. P ff 0.05, 
1.6 DEmiRNAs 靶 基 因 预 测 及 分 析 

利用 TargetFinder 软件 对 DEmiRNAs 进行 靶 基 
因 预 测 ( Allen et al., 2005) ,并 利用 Blast 软件 将 预 
测 丢 基因 序列 与 GO( Gene Ontology ) 数据 库 进行 比 
对 ,获得 靶 基 因 的 注释 信息 ,将 miRNAs 的 节 基 因 按 
照 细胞 组 分 、 分 子 功 能 以 及 生物 学 进程 进行 分 类 。 
采用 基 迪 奥 在 线 分 析 软 件 ( http: // www. omicshare. 
com/tools/ Home/ Index/index. html) 对 DEmiRNAs 的 
靶 基 因 进 行 KEGG pathway 富 集 分 析 。 最 后 ,根据 
DEmiRNA 与 mRNAs 的 靶 向 结合 关系 , 利用 
Cytoscape 软件 构建 miRNAs-mRNAs 的 调控 网 络 。 
1.7 DEmiRNAs 的 茎 环 实时 荧光 定量 PCR ( stem- 
loop qRT-PCR ) 检测 

根据 DEmiRNAs 的 序列 信息 ,利用 DNAMAN 软 





























件 设计 stem-loop 引物 、 正 向 引物 和 通用 反 向 引物 ， 
并 委托 上 海 生 工 生物 工程 公司 合成 相关 引物 。 本 研 
究 使 用 引物 信息 详 见 表 1。 利 用 RNA 抽 提 试剂 盒 
(Axygen 公司 ,美国 ) 分 别提 取 AmCK 和 AmT 的 总 
RNA ,利用 stem-loop 引物 进行 反 转 录 获 得 相应 的 
cDNA 模板 。 常 规 PCR 体系 为 20 uL: PCR mix 
(TaKaRa 公司 ,日 本 ) 10 uL, 无 菌 水 7 uL, E FUR 
引物 (1.67 jumol/L) 及 cDNA 模板 各 1 uL. RMF 
序 : 95*C 3 min; 95*C 30 s, 49% 30 s, 72C 1 min, 
34 个 循环 ; 72% 5 min, PCR 的 产物 经 1. 896 的 琼 
和 糖 凝 胶 电泳 检测 。qPCR 检测 采用 SYBR Green 
法 在 ABI QuantStudio 3 KJEE Œ PCR 系统 (ABI 公 
n] ,美国 ) 中 进行 ,反应 体系 为 20 pL: SYBR Green 
Dye 10 pkL， 上 下 游 引物 (1.67 pmol/L) 以 及 cDNA 
模板 各 1 uL, Rox 0.44 uL, DEPC 水 补 至 20 pL。 
反应 条 件 : 95% 预 变性 1 min; 95% 变性 15 s, 49*C 
延伸 30 s, 40 个 循环 ; 最 后 72% 延伸 45 s。 以 UG 
作为 内 参 , 所 选 mi RNAs 的 相对 表达 量 采 用 2 A 
式 进 行 计 算 。 对 每 个 生物 学 重复 样品 进行 3 次 技术 
重复 。 最 后 通过 Graph Prism 5 软件 对 DEmiRNA 的 
qPCR 结果 和 sRNA-seq 结果 进行 绘图 。 









































X1 本 研究 使 用 的 引物 
Table 1 Primers used in this study 


引物 Primers 


引物 序列 Primer sequences (5' -37) 





AmLOOP-miR-1672-x 
AmLOOP-miR-4968-y 
AmLOOP-miR-251-x 
AmLOOP-miR-9277-y 
AmF-miR-1672-x 
AmF-miR-4968-y 
AmF-miR-251-x 
AmF-miR-9277-y 
AmR 
U6-F 
U6-R 


2 结果 


2.1 意大利 蜜蜂 6 日 龄 幼虫 肠 道 样品 测序 数据 质 
控 与 评估 

本 研究 中 ,AmCK 和 AmT 两 组 幼虫 肠 道 样品 的 
sRNA-seq 分 别 产生 11 171 844 和 13 160 362 条 raw 
reads ,经 严格 过 滤 后 得 到 的 clean tags 数 分 别 为 





CTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGTCCCTTCC 
CTCAACTGGTGTCGTGGAGTCGGCA ATTCAGTTGAGTGCTGCTG 
CTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGTTTCCTGC 
CTCAACTGGTGTCGTGGAGTCGGCAATTCAGTTGAGAGTCGGTA 


CGGTCAGGCCC 
CAGCAGCAGCAG 
ATAAGTACTAGTGCC 
AGGTTCGAGTCC 
CTCAACTGGTGTCGTGGA 
GTTAGGCTTTGACGATTTCG 
GGCATTTCTCCACCAGGCTA 


9 230 496 (84.56% ) 和 10 823 667 (84. 7796) ( x 
2), AmCK 和 AmT 组 内 各 生物 学 重复 之 间 的 
Pearson 相关 性 系数 均 在 0.993 以 上 ,说 明 各 样品 的 
复 性 良好 (图 1) 。 上 述 结果 说 明 本 研究 的 测序 数 
据 质量 良好 ,可 用 于 进一步 分 析 。 
2.2 意大利 蜜蜂 6 日 龄 幼虫 肠 道 样品 DEmiRNAs 
的 预测 及 分 析 

利用 miRDeep2 软件 共 预 测 出 740 个 AmCK vs 























1 期 BO RT: 意大利 蜜蜂 幼虫 肠 道 在 球赛 菌 胁迫 后 
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AmT 的 共有 miRNAs, 其 中 包括 ame-let-7 ( TPM = 
5 152, log;,FC = - 0. 18) 和 ame-bantam ( TPM = 
37 220, log; FC = -0.27) YE ÉY 305 个 mi RNAs 在 
AmT 中 的 表达 量 均 为 下 调 。 进 一 步 筛 选 发 现 AmCK 
vs AmT 中 包含 21 个 DEmiRNAs ,其 中 上 调和 下 调 
miRNAs 的 数量 分 别 为 15 和 6 个 ( 表 3)。 绪 果 显 
示 ,miR-8440-y，miR-3245-y 和 miR-6941-x 上 调 表 
达 幅 度 较 大 ,miR-106-y,，miR-429-y 和 novel-m0003- 
3p 下 调幅 度 较 大 ;此 外 ,DEmiRNAs 中 仅 包含 2 个 
novel miRNAs, 分 别 为 novel-m0018-5p 和 novel- 
m0003-3p ,其 余 19 个 均 为 已 知 miRNAs。 


表 2 意大利 蜜蜂 6 日 龄 幼虫 肠 道 样品 
SRNA-seq 数据 概况 
Table 2 Overview of sRNA-seq datasets of gut 
samples from the 6-day-old larvae of 


Apis mellifera ligustica 








样品 原始 读 段 有 效 标 签 序列 数 ( 占 比 ) 
Samples Number of raw reads Number of clean tags ( Proportion) 
AmCK-1 11 470 909 9 324 373 (82.80% ) 
AmCK-2 10 773 616 9 038 930 (86.63% ) 
AmCK-3 11 271 007 9 328 185 (84.26% ) 
AmT-1 12 974 974 10 718 936 (84.86% ) 
AmT-2 16 781 118 13 672 568 (83.7796 ) 
AmT-3 9 724 993 8 079 498 (85.67% ) 


AmCK; 正常 意大利 蜜蜂 对 照 组 Normal A. m. ligustica as the control 
group; AmT: 球 吉 菌 侵 染 的 意大利 蜜蜂 处 理 组 A. apis challenged A. 


m. ligustica as the treatment group. F |E] The same below. 
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图 1 意大利 蜜蜂 6 日 龄 幼虫 肠 道 样品 的 不 同 
生物 学 重复 间 的 Pearson 相关 性 


Pearson correlations between different biological 








Fig. 1 
repeats within each gut sample from the 6-day-old 


larvae of Apis mellifera ligustica 


表 3 正常 及 球 圳 菌 侵 染 的 意大利 蜜蜂 6 E S 
肠 道 样品 中 差异 表达 miRNAs ( DEmiRNAs) 
信息 统计 
Table 3 Summary of information of differentially 
expressed miRNAs ( DEmiRNAs) in the gut of 
the 6-day-old larvae between the normal and 
Ascosphaera apis challenged 
Apis mellifera ligustica 




















METTA 
DEmiRNAs ID AmCK AmT log, FC Number of 
target genes 

miR-8440-y 0.01 9.92 9.95 2 265 
miR-3245-y 0.01 5.33 9.06 85 
miR-6941 -x 0.01 3.19 8.32 79 
miR-1344-x 0.31 17:53 5.81 431 
miR-6497-x 1.04 9.27 3.16 52 
miR-5119-y 4.57 25:32 2.47 701 
miR-250-x 2.98 15.27 2.36 212 
miR-511-y 1.83 8.35 2.19 148 
ame-miR-6046 1.26 5.47 2.12 497 
miR-8117-y 13.56 52. 82 1.96 704 
miR-3793-x 7.09 26.24 1.89 398 
miR-342-y 14.60 48.33 1.73 249 
miR-4448-y 265.71 576.57 1.12 88 
novel-m0018-5p 2.06 4.37 1.09 72] 
ame-miR-3759 174.09 350.19 1.01 221 
miR-81-y 25.45 8.80 -1.53 581 
miR-1-z 15.27 0.80 -4.26 645 
ame-miR-927b 1.04 0.01 -6.70 1 922 
novel-m0003 -3p 1.04 0.01 -6.70 297 
miR-429-y 2.08 0.01 -7.70 993 
miR-106-y 4.42 0.01 -8.79 506 


log; FC; 以 2 HJE miRNA 的 相对 变化 倍数 的 对 数值 log; Fold change. 


2.3 ”意大利 蜜蜂 6 日 龄 幼虫 肠 道 样 品 DEmiRNAs 
的 靶 基 因 预 测 及 功能 注释 

AmCK vs AmT 中 DEmiRNAs 共 预 测 出 6 759 个 
靶 基 因 , 其 中 上 调和 下 调 miRNAs 可 分 别 彼 向 结合 
4 632 和 3 909 EIEE GO 数据 库 注释 结果 显 
示 , 上 调 miRNAs 的 3 503 个 靶 基 因 共 涉及 40 个 GO 
terms , 富 集 靶 基 因数 最 多 的 分 别 是 结合 (1 016 个 )、 
细胞 进程 (941 个 ) .代谢 进程 (762 个 ) . 单 组 织 进程 
(695 ^) 催化 活性 (605 个 )、 细 胞 膜 (424 个 ) 、 细 
胞 膜 组 件 (419 个 ) 细胞 (391 个 ) 、 细 胞 组 件 (391 
个 ) 定位 (285 个 ) 等 (图 2: A); 下 调 miRNAs 的 
3 252 个 靶 基 因 可 注释 到 39 个 GO terms, 富 集 靶 基 
因数 最 多 的 分 别 是 结合 (870 个 ) 细胞 进程 (828 
个 ) 、 单 组 织 进程 (727 个 ) 代谢 进程 (718 个 ) 催化 
活性 (600 个 ) 细胞 (422 个 ) 细胞 组 件 (422 个 )、 
细胞 膜 (357 个 ) 、 细 胞 膜 组 件 (348 个 ) 生物 学 调控 
(312 个 ) 等 (图 2: B), 
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结构 分 子 活性 Structural molecule activity 

转录 因子 活性 ， 蛋 白 结合 Transcription factor activity, protein binding 
抗 氧化 活性 Antioxidant activity 
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细胞 进程 Cellular process 

代谢 进程 Metabolic process 

单 组 织 进程 Single-organism process 

定位 Localization 

生物 学 调控 Biological regulation 

应 激 反应 Response to stimulus 

信号 Signaling 

细胞 成 分 组 织 或 生物 合成 Cellular component organization or biogenesis 
多 细胞 生物 进程 Multicellular organismal process 
生物 粘 附 Biological adhesion 

发 育 进程 Developmental process 

行为 Behavior 

免疫 系统 进程 Immune system process 
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分 子 传感器 活性 Molecular transducer activity 
转运 子 活性 Transporter activity 

信号 传感器 活性 Signal transducer activity 

子 活 性 Nucleic acid binding transcription factor activity 
分 子 功能 调节 器 Molecular function regulator 
结构 分 子 活性 Structural molecule activity 

抗 氧化 活性 Antioxidant activity 

子 活 性 ， 蛋 白 结合 Transcription factor activity, protein binding 
电子 载体 活性 Electron carrier activity 
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图 2 EXKI 








侵 染 的 意大利 蜜蜂 6 日 龄 幼虫 肠 道 样品 中 差异 表达 mi RNAs ( DEmiRNAs) B SEJEDS GO 分 类 


Fig. 2 GO categorization of the target genes of differentially expressed miRNAs ( DEmiRNAs) in the gut of the 6-day-old 
larvae between the normal and Ascosphaera apis challenged Apis mellifera ligustica 
A: EJ] miRNAs 的 靶 基 因 Target genes of up-regulated miRNAs; B; 下 调 miRNAs 的 对 基因 Target genes of down-regulated miRNAs. 


进一步 对 DEmiRNAs fj 58 3t pj Jt f; KEGG 
pathway 富 集 分 析 ,结果 显示 上 调 miRNAs 的 车 基 因 


?一口 





可 注释 到 104 条 代谢 通路 ,其 中 富 集 数 最 多 的 是 
Wnt 信号 通路 (57 个 ) 内 吞 作用 (56 A) Hippo fri 


























郭 SRI 意大利 蜜蜂 幼虫 肠 道 在 球 吉 菌 胁迫 后 期 的 差异 
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号 通路 (43 个 ) 磷脂 酰 肌 醇 信号 系统 (43 个 )、 味 叭 
代谢 (35 AP) 、 神 经 活性 的 配 体 - 受 体 相 互 作用 (33 
^r) 磷酸 肌 醇 代 谢 (32 个 )、Fox0 信号 通路 (31 
A) Hedgehog 信号 通路 (23 个 ) 和 层 夜 节律 (22 
个 ) 等 ( 表 4)。 

下 调 miRNAs 的 靶 基 因 注 释 到 99 条 代谢 通路 ， 




































































包括 Wnt 信号 通路 (81 个 ) 、 光 传导 (50 A) Hippo 信 
号 通路 (38 个 ) mRNA 监视 通路 (35 个 ) 泛 素 介 导 的 
蛋白 水 解 (35 个 )、 碳 代谢 (32 个 ) 内 质 网 内 蛋白质 
加 工 (30 个 ) , FoxO 信号 通路 (27 个 ) A fe EH C25 
个 )、 味 叭 代谢 (23 个 ) 等 ( 表 5)。 此 外 ,分 别 有 8 和 
6 个 对 基因 富 集 在 Jak-STAT 和 MAPK 信号 通路 。 







































































mA 正常 及 球 圳 菌 侵 染 的 意大利 蜜蜂 6 日 龄 幼虫 肠 道 样品 中 上 调 miRNAs 的 靶 基 因 富 集 数 前 20 个 代谢 通路 
Table 4 Top 20 metabolic pathways enriched by the target genes of the up-regulated miRNAs in the gut from 
the 6-day-old larvae of the normal and Ascosphaera apis challenged Apis mellifera ligustica 















































代谢 通路 靶 基 因 富 集 数 代谢 通路 ID 
Metabolic pathways Number of enriched target genes Pathway ID 
Wnt 信号 通路 Wnt signaling pathway 57 ko04310 
WAE Endocytosis 56 ko04144 
Hippo 信和 号 通路 Hippo signaling pathway 43 ko04391 
磷脂 酰 肌 醇 信号 系统 Phosphatidylinositol signaling system 43 ko04070 
味 叭 代谢 Purine metabolism 35 ko00230 
神经 活性 的 配 体 - 受 体 相 互 作 用 Neuroactive ligand-receptor interaction 33 ko04080 
磷酸 肌 醇 代谢 Inositol phosphate metabolism 32 ko00562 
FoxO 信和 号 通路 FoxO signaling pathway 31 ko04068 
Hedgehog 信和 号 通路 Hedgehog signaling pathway 23 ko04340 
层 夜 节律 Circadian rhythm 22 ko04711 
mTOR 信号 通路 mTOR signaling pathway 21 ko04150 
溶 酶 体 Lysosome 20 ko04142 
光 传 导 Phototransduction 19 ko04745 
淀粉 与 荐 糖 代谢 Starch and sucrose metabolism 19 ko00500 
泛 素 介 导 的 蛋白 水 解 Ubiquitin-mediated proteolysis 19 ko04120 
mRNA 监视 通路 mRNA surveillance pathway 18 ko03015 
氨基 糖 与 核 苷 酸 糖 代谢 Amino sugar and nucleotide sugar metabolism 17 ko00520 
Notch 信和 号 通路 Notch signaling pathway 17 ko04330 
内 质 网 内 和 蛋白质 加 工 Protein processing in endoplasmic reticulum 17 ko04141 
碳 代 谢 Carbon metabolism 16 ko01200 


表 5 


正常 及 球 圳 菌 侵 染 的 意大利 富 蜂 6 日 龄 幼虫 肠 道 样品 中 下 调 miRNAs 的 靶 基 因 富 集 数 前 20 个 代谢 通路 


Table 5 Top 20 metabolic pathways enriched by the target genes of the down-regulated miRNAs in the gut from 
the 6-day-old larvae of the normal and Ascosphaera apis challenged Apis mellifera ligustica 


代谢 通路 


Metabolic pathways 











Wnt 信和 号 通路 Wnt signaling pathway 

光 传 导 Phototransduction 

Hippo 信号 通路 Hippo signaling pathway 
mRNA 监视 通路 mRNA surveillance pathway 
泛 素 介 导 的 和 蛋白 水 解 Ubiquitin-mediated proteolysis 

碳 代 谢 Carbon metabolismn 

内 质 网 内 和 蛋白质 加 工 Protein processing in endoplasmic reticulum 
FoxO 信和 号 通路 FoxO signaling pathway 

WFE 
味 叭 代谢 Purine metabolism 

组 氮 酸 代谢 Histidine metabolism 

磷酸 肌 醇 代谢 Inositol phosphate metabolism 
剪接 体 Spliceosome 

'TGF-beta f, 
磷脂 酰 肌 本 
edgehog 
基 糖 与 核 苷 酸 糖 代谢 Amino sugar and nucleotide sugar metabolism 
虱 轴 的 形成 Dorsoventral axis formation 

Wi Phagosome 


神经 活性 的 配 体 - 受 体 相 互 作 用 Neuroactive ligand-receptor interaction 
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靶 基 因 富 集 数 代谢 通路 ID 
Number of enriched target genes Pathway ID 
81 ko04310 
50 ko04745 
38 ko04391 
35 ko03015 
35 ko04120 
32 ko01200 
30 ko04141 
27 ko04068 
25 ko04144 
23 ko00230 
22 ko00340 
21 ko00562 
21 ko03040 
20 ko04350 
19 ko04070 
18 ko04340 
17 ko00520 
17 ko04320 
17 ko04145 
16 ko04080 
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2.4 ”意大利 蜜蜂 6 日 龄 幼虫 肠 道 样品 免疫 相关 个 与 细胞 免疫 中 泛 素 介 导 的 蛋白 水 解 相 关 的 


DEmiRNAs 的 调控 网 络 构建 及 分 析 

对 富 集 在 丝氨酸 蛋白 酶 和 泛 素 介 导 的 蛋白 水 解 
上 的 靶 基 因 及 与 其 结合 的 miRNAs 进行 可 视 化 , 结 
果 显 示 17 个 DEmiRNAs 靶 向 结合 116 个 与 体液 免 
疫 中 丝 氮 酸 和 蛋白 酶 相关 的 mRNAs (图 3: A), 
miRNAs 与 mRNAs 之 间 形 成 1 个 较 大 和 2 个 较 小 的 
复杂 调控 网 络 , 其 中 结合 mRNAs 数 最 多 的 
DEmiRNAs 为 ame-miR-927b, miR-8840-y 和 miR- 























mRNAs 被 14 个 DEmiRNAs 靶 向 结合 (图 3: B) , 相 
互 之 间 形 成 1 个 较 大 的 和 5 个 较 小 的 调控 网 络 ,其 
中 ame-miR-927b , miR-429-y 和 miR-8440-y 靶 向 结合 
的 mRNAs 最 多 ,分 别 达 到 21, 13 和 13 个。 此 外 ， 
miR-5119-y, mi R-3793-x , mi R-429-y 和 ame-miR-927b 
可 分 别 靶 向 结合 10, 10, 6, 2 0E 4 Te 2RUIBL US] T 388 
路 上 的 mRNAs( XM, 006565561.2, XM. 016914484. 1, 
XM. 393703.6, XM. 016911719. 1 和 XM, 016911720. 1 


8117-y ,分 别 结合 46, 43 4132 个 mRNAs。 共 有 54 等 ) 。 
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图 3 





富 集 在 意大利 蜜蜂 幼虫 肠 道 免疫 相关 通路 上 的 差异 表达 miRNAs ( DEmiRNAs) 的 调控 网 络 


Fig. 3 Regulation networks of differentially expressed miRNAs ( DEmiRNAs) enriched in immune-related pathways in 
the larval gut of Apis mellifera ligustica 
A: 调控 丝氨酸 蛋白 酶 的 DEmiRNAs 的 调控 网 络 Regulation networks of DEmiRNAs regulating serine protease; B; 调控 泛 素 介 导 的 蛋白 水 解 的 


DEmiRNAs 的 调控 网 络 Regulation networks of DEmiRNAs regulating ubiquitin-mediated proteolysis. 绿色 圆 形 代 表 mRNAs, 红色 三 


Green circles indicate mRNAs, while red triangles indicate miRNAs. 


进一步 分 析 发 现 , ame-miR-927b, miR-5119-y 
和 miR-429-y 共同 参与 调控 上 述 3 种 免疫 相关 通 
路 ;而 miR-3793-x 同时 参与 调控 细胞 凋 亡 和 泛 素 介 
导 的 蛋白 水 解 ;ame-miR-3759 , ame-miR-6046, miR- 











代表 miRNAs。 














106-y, miR-1334-x, miR-250-x, miR-8117-y, miR- 
81-y, miR-8440-y 和 novel-m0018-5p 在 内 的 9 个 
miRNAs 同时 参与 调控 泛 素 介 导 的 蛋白 水 解 和 丝 氮 
ME A M; mi R-6497 x 仅 在 泛 素 介 导 的 蛋白 水 解 中 
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发 挥 调 控 作 用 , miR-1-z，miR-3245-y，miR-342-y， 
miR-511-y 和 novel-m0003-3p Atl ASI ZARE EH 
酶 的 调控 。 
2.5 意大利 蜜蜂 6 日 龄 幼虫 肠 道 样品 DEmiRNAs 
的 stem-loop RT-qPCR 的 验证 

为 整体 验证 测序 数据 的 准确 性 ,从 所 有 显著 差 
异 和 非 显著 差异 的 740 个 miRNAs 中 随机 挑选 2 个 
上 调 miRNAs ( miR-251-x 和 miR-9277-y ) 和 2 个 下 
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图 4 正常 及 球 吉 菌 侵 染 的 意大利 蜜蜂 














调 miRNAs ( miR-1672-x 和 miR-4968-y ) 进行 stem- 
loop RT-PCR 验证 ,电泳 结果 显示 上 述 4 个 
DEmiRNAs 在 AmCK 和 AmT 中 均 可 扩 增 出 符合 预期 
的 目的 片段 (图 4: A, B, ,下 )。 进 一 步 的 RT-qPCR 
结果 显示 ,miR-251-x，miR-1672-x，miR-9277-y 和 
miR-4968-y 的 表达 量变 化 趋势 与 测序 数据 中 相应 的 
DEmiRNAs 的 表达 量变 化 趋势 一 致 (图 4: C, D, G, 
H) 。 上 述 结果 说 明 本 研究 的 测序 结果 真实 可 靠 。 
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6 日 龄 幼虫 肠 道 样品 中 差异 表达 miRNAs ( DEmiRNAs) 的 
stem-loop PCR 和 RT-qPCR 验证 





Fig. 4 Stem-loop PCR and RT-qPCR validation of differentially expressed miRNAs ( DEmiRNAs) in the gut of 
the 6-day-old larvae between the normal and Ascosphaera apis challenged Apis mellifera ligustica 
A, B, E, F; DEmiRNAs 的 PCR 产物 的 琼脂 糖 凝 胶 电 泳 Agarose gel eletrophoresis of the PCR products of DEmiRNAs; C, D, G, H; DEmiRNAs 的 


qPCR 结果 qPCR results of DEmiRNAs. 





3 讨论 





昆虫 的 miRNAs 可 被 病原 菌 诱导 或 抑制 ,从 而 呈现 
上 调 或 下 调 表达 ( 赵 新 慧 , 2013) 。 相 比 于 果 蝇 和 家 
大 等 模式 昆虫 ,蜜蜂 的 miRNA 研究 较为 滞后 ,并 主 
































MIRNA 作为 转录 后 水 平 的 关键 调控 因子 ,在 昆 
虫 的 免疫 防御 过 程 中 发 挥 重 要 作用 ( 吴 萍 , 2016), 














要 集中 在 生长 发 育 (Ashby et al., 2016) 和 级 型 分 化 
(BRIE, 2012) 等 方面 ,关于 miRNA 在 蜜蜂 免疫 防御 
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过 程 中 的 作用 研究 进展 缓慢 。Lourenco 等 (2013 ) 发 
现 被 芯 黄 微 球菌 Micrococcus luteus 和 粘 质 沙 雷 氏 菌 
Serratia marcescens 感染 后 3 西方 蜜蜂 工蜂 的 ame-let- 
7，ame-miR-13a，ame-miR-2 和 ame-miR-92a 表达 量 
上 调 , 而 ame-bantam, ame-miR-12 和 ame-miR-184 
等 8 个 miRNAs 表达 量 下 调 ,作者 推测 DEmiRNAs 
在 蜜蜂 免疫 系统 中 发 挥 潜在 的 调控 作用 。 本 研究 利 
用 sRNA-Seq 技术 对 正常 及 球 早 菌 胁迫 的 意 蜂 6 日 
龄 幼虫 肠 道 进 行 测序 ,预测 出 15 个 上 调 miRNAs 和 
6 个 下 调 miRNAs ,表明 球 才 菌 胁迫 可 引起 意 蜂 幼虫 
肠 道 miRNAs 的 变化 ,推测 DEmiRNAs 可 参与 到 宿 
主 的 胁迫 应 答 过 程 。 

本 研究 中 ,对 于 上 调 miRNAs 的 靶 基 因 , 富 集 在 
细胞 进程 (941 个 ) .代谢 进程 (762 个 ) 细胞 膜 (424 
个 ) 等 GO terms 上 的 数量 , 均 高 于 富 集 在 上 述 terms 
的 下 调 miRNAs 的 靶 基 因数 ,表明 意 蜂 幼虫 肠 道 的 
新 陈 代谢 ,细胞 生命 活动 在 球 守 菌 胁迫 后 期 均 受 到 
一 定 程 度 的 抑制 。Wnt 信和 号 通路 和 FoxO 信和 号 通路 
可 参与 到 哺乳 动物 生殖 系统 的 形成 与 发 育 (Zheng et 
al., 2010; Edmonds et al., 2010) , E Eb E Hj fi A LA 
及 寿命 调节 等 方面 也 发 挥 重 要 的 调控 作用 ( 陈 伟 和 
徐 卫 华 , 2015; 刘 艳 玲 ,2016 ) 。Hippo 信号 通路 和 
Notch 信号 通路 的 共同 作用 可 影响 果 蝇 的 器 官 发 育 
( Barry and Camargo, 2013) , Hedgehog 信和 号 通路 可 
参与 到 B. mori BG 28 8L ^3 16 zr B9 e s (3E 36, 
2014)。 本 研究 发 现 , 在 Wnt, Hippo, FoxO, 
Hedgehog 以 及 Notch 信号 通路 上 分 别 富 集 了 57, 
43,31, 23 和 17 个 上 调 miRNAs WEAK 4) ， 
表明 这 些 上 调 miRNAs 广泛 参与 了 上 述 信号 通路 的 
负 调 探 过 程 , 意 蜂 幼虫 肠 道 吉 官 生长 .组 织 发 育 相关 
的 信号 通路 受到 球 喜 菌 的 抑制 。 此 外 还 发 现 ,下 调 
miRNAs 同样 参与 调控 意 蜂 幼虫 肠 道 的 生长 发 育 和 
新 陈 代谢 过 程 如 碳 代 谢 (32 A) UURS ERE FORE C21 
个 ) 和 青 腹 轴 的 形成 (17 个 ) 等。 本 研究 中 , 分 别 有 
8 和 6 个 下 调 miRNAs 的 靶 基 因 富 集 在 Jak-STAT 和 
MAPK 等 体液 免疫 相关 的 信号 通路 ,35 和 25 个 下 
调 miRNAs 的 靶 基 因 注 释 到 泛 素 介 导 的 蛋白 水 解 和 
内 吞 作 用 等 细胞 免疫 相关 的 通路 , 推测 这 些 
miRNAs 通过 下 调 表达 量 ,部 分 解除 对 相关 免疫 基 
的 抑制 ,从 而 增强 对 球 圳 菌 的 抵御 。 

昆虫 肠 道 可 通过 产生 抗菌 肽 和 活性 氧 抵御 病原 
微生物 的 侵 染 ( 陈 康 康 , 2016) 。 丝 氨 酸 和 蛋白酶 参与 
了 昆虫 体液 免疫 过 程 中 Toll 免疫 信号 通路 和 胞 内 酚 
氧化 酶 原 的 激活 ,在 血 淋 巴 黑 化 和 抗菌 肽 合成 等 免 



















































































































































































































































































疫 防 御 过 程 中 起 关键 作用 ( 唐 琦 等 , 2017) 。 昆 虫 体 
内 含有 病原 识别 的 受 体 蛋 白 , 可 识别 真菌 来 源 的 B- 
1,3- 和 区 聚 糖 ,从 而 激活 细胞 内 丝氨酸 蛋白 酶 的 活性 
以 应 对 真菌 的 侵 染 ( Amparyup et al., 2013) ,例如 亚 
JI SE DR. Ostrinia furnacalis 体内 的 丝氨酸 和 蛋白酶 
SP1 和 SPI3 可 介 导 黑 化 反应 以 应 对 真 戎 的 入 侵 
(Chu et al.，2015)。 本 研究 中 ,共有 17 个 
DEmiRNAs 在 意 蜂 幼虫 肠 道 响应 球 吉 戎 胁迫 过 程 中 
参与 调控 丝氨酸 和 蛋白酶 通路 ,其 中 11 个 上 调 
miRNAs( miR-250-x, miR-8117-y 和 miR-8440-y 等 ) 
丢 向 结合 93 个 mRNAs,6 个 下 调 miRNAs( miR-1-z, 
miR-106-y 和 ame-miR-927b 等 ) 靶 向 结合 74 个 
mRNAs; 多 数 miRNAs 上 调 表 达 水 平 以 加 强 对 苇 
mRNAs 的 抑制 作用 ,少数 miRNA 下 调 表 达 水 平 以 
减弱 对 靶 mRNAs HIMI; ERRAK HERE KTE 
迫 后 期 对 宿主 的 丝氨酸 蛋白 酶 通路 产生 影响 ,二 者 
之 间 存 在 复杂 的 互 作 , 但 背后 的 分 子 机 制 需 要 进 
步 探 明 。 

miRNAs 具有 高 度 保守 性 、 组 织 特异 性 和 表达 
时 序 性 (He and Hannon, 2004), Garbuzov 等 
(2010) 发 现 let-7 可 抑制 果 晶 体内 抗菌 肽 的 翻译 过 
程 ,从 而 调节 其 天 然 免 疫 反应 。Etebari 等 (2013 ) 发 
现 被 弯 尾 姬 蜂 Diadegma semiclausum 寄生 的 小 菜 蛾 
Plutella xylostella 幼虫 体内 let-7 的 表达 量 明 显 提 高 ， 
并 参与 了 宿主 的 免疫 调控 。 同 源 mi RNAs 在 不 同 昆 
虫 的 免疫 防御 中 可 能 发 挥 类 似 的 调控 作用 。 本 研究 
发 现 , 意 蜂 幼 虫 肠 道 的 ame-let-7 在 球 申 菌 胁迫 后 期 
表达 量 降低 ,推测 宿主 通过 下 调 ame-let-7 的 表达 水 
平 减弱 对 靶 基 因 的 抑制 作用 ,从 而 激活 相关 免疫 通 
路 ,以 应 对 球赛 菌 的 侵 染 。Su 等 (2017 ) 发 现 白 纹 伊 
Ty. Aedes albopictus 在 摄食 含 登 革 病 毒 CDENV ) 85 ifi 
液 后 其 中 肠 内 miR-1 和 miR-4448 表达 量 分 别 呈 上 
调和 下 调 表 达 。 本 研究 发 现 , 意 蜂 幼 虫 肠 道 的 mi R- 
1-z 和 miR-4448-y 的 表达 量 在 球 上 时 阔 胁迫 后 期 分 别 
下 调和 上 调 表 达 , 表 明 miR-1-z 和 miR-4448-y 涉及 
宿主 对 球 蝇 菌 的 胁迫 应 答 , 此 外 ,本 研究 中 miR-1-z 
和 miR-4448-y 与 白 纹 伊 蚊 中 的 同 源 miRNA 表达 趋 
势 情况 相反 ,推测 同 源 mi RNAs 在 不 同 物种 中 可 能 
发 挥 不 同 的 作用 ,但 这 种 差异 产生 的 分 子 机 制 仍 需 
要 进一步 研究 。 细 胞 凋 亡 在 机 体 细胞 免疫 和 抵御 病 
原 微生物 入 侵 的 过 程 中 扮演 重要 角色 ( 包 希 艳 等 ， 
2017)。 家 不 在 被 家 盔 核 型 多 角 体 病毒 ( BmNPV ) 侵 
染 的 早期 ,体内 的 促 凋 亡 基 因 呈 上 调 表 达 , 清除 有 
害 、 宛 余 的 细胞 ， 以 保护 其 他 细胞 免 受 伤害 (Wu et 
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al., 2013) Campbell 等 (2014 ) 研究 发 现 miR-5119- 
5p 参与 了 埃及 伊 蚊 Aedes aegypti 响应 DENV 侵 染 的 
表达 调控 。 本 研究 中 , 意 蜂 幼 虫 肠 道 的 miR-5119-y 
在 球 宫 菌 胁迫 后 期 表达 量 上 调 ,其 靶 向 结合 的 10 个 
HIJE] ( XM _006565561. 2, XM. 016914484. 1 和 
XM_393703.6 等 ) 富 集 在 细胞 凋 亡 ,推测 球 沾 菌 在 胁 
迫 后 期 能 够 通过 影响 宿主 miR-5119-y 的 表达 ,在 一 
定 程 度 上 抑制 细胞 凋 亡 的 产生 ,从 而 促进 自身 在 宿 
主体 内 的 增殖 。 

miR-927 作为 一 种 高 度 保守 的 miRNA ,参与 到 
意 蜂 幼虫 至 成 年 蜜蜂 羽化 出 房 过 程 中 赔 皮 激素 受 体 
蛋白 的 调控 ,在 相关 受 体 被 敲 除 的 意 蜂 体内 ,ame- 
miR-927a-5p 和 ame-miR -927b-5p 均 下 调 表达 (Mello 
et al., 2014) 。 此 外 ,在 被 DENV 感染 的 C6/36 细胞 
内 ,miR-927 显著 上 调 表 达 , 其 靶 基 因 可 能 参与 了 
Jak-STAT 级 联 反应 , Toll 信号 通路 以 及 相关 蛋白酶 、 
ATP 酶 的 调控 ( Avila-Bonilla et al., 2017)。 本 研究 
中 , 意 蜂 幼虫 ame-miR-927b 的 表达 量 在 球 圳 菌 侵 染 
后 期 , 呈 显 著 下 调 (log,FC = -6.7) 并 在 与 丝氨酸 蛋 
白 酶 泛 素 介 导 的 蛋白 水 解 相关 的 调控 网 络 中 , 靶 向 
结合 mRNAs 数量 均 为 最 多 (46 和 21 个 )。 上 述 结 
果 表 明 ame-miR-927b 在 意 蜂 幼虫 肠 道 响应 球 直 菌 
胁迫 的 免疫 防御 中 具有 关键 的 调控 功能 。 下 一 步 可 
通过 合成 相应 的 miRNA 类 似 物 (mimic) 和 miRNA 
抑制 物 (inhibitor) 对 其 进行 过 表达 和 敲 减 ,从 而 深入 
研究 其 功能 。 本 研究 仅 对 参与 调控 丝氨酸 蛋白 酶 、 
泛 素 介 导 的 蛋白 水 解 的 DEmiRNAs 及 其 靶 向 
mRNAs 构建 调控 网 络 ,后 续 研 究 将 对 涉及 调控 其 他 
免疫 通路 或 代谢 通路 的 miRNAs 及 其 靶 向 mRNAs 
构建 调控 网 络 ,进一步 第 选 居 于 核心 位 置 的 关键 
miRNAs, 

Zk EIIE , BIFX AE 76 A ERIE E A E] E E 
6 日 龄 幼虫 肠 道 样品 进行 高 通 量 测序 ,对 宿主 的 
miRNAs 表达 谱 和 DEmiRNAs 进行 深入 分 析 。 研 究 
结果 不 仅 为 意 蜂 幼虫 肠 道 响应 球面 菌 胁迫 的 关键 调 
控 miRNAs 的 筛选 及 功能 研究 提供 了 重要 的 信息 和 
线索 ,也 为 意 蜂 幼虫 - 球 圳 菌 互 作 机 制 的 前 明 打 下 了 
基础 。 
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